Genuine three-partite entangled states with a local hidden variable model 
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We present a family of three-qubit quantum states with a bàsic local hidden variable model. Any 
von Neumann measurement can be described by a local model for these states. We show that some 
of these states are genuine three-partite entangled and also distillable. The generalization for larger 
dimensions or higher number of parties is also discussed. As a byproduct, we present symmetric 
extensions of two-qubit Werner states. 
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I. INTRODUCTION 

One of the most striking characteristics of quantum 
mechanics is nonlocality. If quantum mechanics could 
be described by a local hidden variable model (LHV) 
then the vàlues measured for multi-particle correlations 
could be reproduccd assuming that all measurable single- 
particle operators had already a value before the mea- 
surement. Bell showed that there are quantum states for 
which the many-body correlations cannot be explained 
based on this assumption However, such correlations 
arise only for some entangled quantum states while for 
separable states the correlations can always be mimicked 
by a LHV model 0. 

Proving that the measurement results on a quantum 
state cannot be obtained from a LHV model is done by 
finding a Bell incquality which is violated by the state 
[H. However, this is difficult since the determination of 
all Bell incqualities is a computationally hard problem 
3]. To prové that any measurement on a given state 
can be described by a LHV model is perhaps even more 
challenging. This is because in order to do that one has to 
find a LHV model for any number of arbitrary operators 
measured at each party. 

Due to the difficulty of the problem, LHV models have 
a quite limited literature. The first and most fundamen- 
tal result of the subject was presented by Werner in Ref. 
Q. He described a LHV model for arbitrary von Neu- 
mann measurements for some U<S>U symmetric bipartite 
states . For the qubit case these are of the form 



Pw =p\ip I + (1 -p) 



4' 



(1) 



where \ip~) = (|01) - |10))/v^ is the singlet state. In the 
same paper Werner also gave the first modern definition 
of quantum entanglement, thus distinguishing it from 
nonlocality. Indeed, states pw are entangled for p > 1/3 
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and local for p < 1/2 2]. It is hard to overestimate 
the importance of these results for the development of 
quantum information science. Later, Barrett obtained a 
model for general measurements, also called positive op- 
erator valued measures (POVMs), for a subset of Werner 
states Q. LHV models were also constructed for finite 
number of settings for states with positive partial trans- 
pose exploiting symmetric extensions Q. Apart from 
their fundamental interest, these results are also relevant 
from a quantum information theory viewpoint. Simulat- 
ing entanglement by classical means (e.g., Ref. sheds 
light on the power of entanglement as information re- 
source. In this context, those quantum states for which 
the correlations can be reproduced by a LHV model are 
useless for communication tasks, since they do not pro- 
vide any advantage over shared classical randomness 0] • 

New and interesting open qüestions on the relation- 
ship between nonlocality and entanglement appear in the 
multipartite scenario. Recali that multipartite entangle- 
ment is known to be inequivalent to bipartite entangle- 
ment 8]. Moreover, genuine multipartite entanglement 
is the property most often detected in experiments (e.g., 
Ref. Q). We know that the Bell inequality violation rc- 
quired for genuine multipartite entanglement (i.e., when 
all parties are entangled with each other Í10| ) increases 
exponentially with the number of parties Hence 
one could expect that entanglement of this type for large 
enough number of parties provides a sufncient condition 
for a state to be nonlocal. However, beyond the bipartite 
case, the connection between nonlocality and entangle- 
ment remains largely unexplored. Indeed, LHV models 
for multipartite entangled systems are still missing. 

In this paper, we present a one-parameter family of 
three-qubit states whose correlations for von Neumann 
measurements can be reproduced by a LHV model. Thus 
these states do not violate any Bell inequality. Then, we 
prové that, remar kably, some of these states have genuine 
three-qubit entanglement [ï(| and we also show that they 
are distillable. To our knowledge, these are the first ex- 
amples of genuine multipartite entangled states allowing 
for a local description. The generalization of the con- 
struction to other situations, more parties or higher di- 
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FIG. 1: Schematic representation of our two-qubit hidden 
variable model. Party B receives a Standard Bloch vector. 
The dashed arrow on the left hand side points opposite to 
this vector. Party A receives a vector pointing to one of the 
eight vertices of a cube tangent to the Bloch sphere. The 
vèrtex is chosen such that the overlap with the dashed vector 
is maximal. 

mensional systems, is also discussed. 

Before proceeding, let us introduce the notation. We 
denote von Neumann measurements on n parties A, B, 
C, ... as Ma, Mb, Mq, ■■■■ The spectral decomposition of 
Ma is given as Ma — J2Í=i a kPk- In the case of qubits, 
that we mostly consider in this work, ot\ — +1 and «2 = 
— 1, while Ma — n A -a, where tia is the normalized vector 
dcfimng the direction of the von Neumann measurement 
and a is the vector of Pauli matrices, er = (a x , a y , a z ). 

II. TWO-QUBIT CASE 

The key point for the construction of our LHV model 
for three qubits is an alternative derivation of Werner's 
original result for two qubits. This new derivation has 
the advantage of being easily generalizable to the case of 
three qubits. Consider the two-qubit operators 

P {2 ' c) := / M{du)g u ® p m (2) 

where 

Qlu ■= \[^- c ^2 SÍgn((o- fe ) w )cr fc ], 

k— x,y,z 

:= \uj)(lj\. (3) 

Here \lo) is a two-element state vector and M is the 
unique probability measure invariant under all single- 
qubit unitary rotations. Direct calculation shows that 
^(2,c) are Werner states ([T]) with p = c/2. Based on this 
construction, the following statement can be made 
Theorem 1 /£/: There exists a LHV model for von Neu- 
mann measurements on states p^ 2 ^ for c < 1. 

Before starting the proof, let us explain the intuition 
behind it. First of all, note that one can restrict the 



analysis to c = 1. In this case, the decomposition ([2]) 
can be understood as a sort of local model for which 
party B receives a Standard Bloch vector, n^, while A 
receives one of the vectors pointing to the vertices of a 
cube. This cube is tangent to the Bloch sphere for the 
Bloch vectors which point to the directions ±x, ±y and 
±z (see Fig.Q]). Actually, A receives the vector with max- 
imum overlap with — n^. Using the Standard trace rule 
(Ma <8 Mb) = J M{duj)Tr(M A e w )Tr(M B p u ), we have 
a LHV model for M A = ±a x / y / z and arbitrary Mg. 
The choice of Ma is restricted since if party A chooses 
other operators to measure, for some ui she would ob- 
tain \(A1a)q^\ > 1- We then say that party B has a 
physical qubit, while A is receiving a non-physical Bloch 
vector. As we have already said, state p^ 2 ' 1 ^ is a Werner 
state, i.e., it is invariant under transformations of the 
form U ®U where U is an arbitrary unitary matrix. Us- 
ing this symmetry, we can construct a LHV model for all 
measurements. The detailed proof goes as follows: 

Proof of Theorem 1. The goal is to find a LHV model 
for the state @ with c = 1, that is to write its correla- 
tions as 

Tr(M A ® M B p (24) ) = I M{du){M A ) u (M B ) u , (4) 

J u 

where (Ma/b)u> are the expectation vàlues of Ma/b if 
the value of the hidden variable is u>, and we require 
\(Ma}uj\,\(Mb)lü\ < 1- Idcntifying the sub-enscmble Ín- 
dex lo in Eq. with the hidden variable in Eq. (QJ , one 
has a LHV model with 

(M A ) U = Tr(M AQa ) = -~Tr[M A ^ sign^^K], 

k— x,y,z 

(M B ) U = Tr{M BPu ) = Tc{Mb\u)(u\). (5) 

It is clear that this model works only if M A — ±a x / y / z . 

Now we modify our LHV model in order to allow ar- 
bitrary operators Ma of the type fÏA&. Such an operator 
can be written in the form 

M A = U\a z U A . (6) 

We can take advantage of the invariance of Werner 
states under transformations of the form U <g> U, so 
(M A ®M B ) = (a z A) (g>M' B ), where M' B = U A M B U A . 
Hence 

(M A ) U = Tr(a z g u ) = -sign[Tr(a z \uj)(oj\)}, 
(M B ) U = Tr(M' B p u ) = Tr(U A M B U] í \u){w\). (7) 

Indeed, substituting Eq. ([7]) into Eq. ((4]) reproduces all 
two-qubit correlations (Ma® Mb). However, Eq. (J7J) is 
not a LHV model for arbitrary Ma and Mb yet. This is 
because (Mb)u depends on Ua, i-e, it depends on what 
operators are measured on party A. This dependence 
can be removed by defming \lo') — u\\oj). We obtain the 
desired LHV model with lo' as the hidden variable given 
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{M A )„> = -sign[Tr(M A |a/)(a/| 
(Mb)^ = TV(MbK)<w'|). 



(8) 



This can trivially be extended to arbitrary Ma UM ■ nc 
can recognize now Werner's model @ where 



1 if (u\P k \uj) < (u\Pi\w) for alH ^ k 



(9) 



otherwise 
for party A while 

(Af B ) w = Tr(M s |w)(w|) (10) 
for B. This finishes the proof of Theorem 1. 

III. THREE-QUBIT CASE 

The previous construction is, in principle, easy to gen- 
eralize to morè parties as follows 



■i n < c ) — 



ec 2 >|=i 



M(du)g u ®p^ n - 1 l (11) 



Here, party A receives again a non-physical Bloch vector, 
while B, C,... get a Standard qubit. We are now in the 
position to prové the main result of this work. 
Theorem 2: There exists a LHV model for von Neumann 
measurements on states p^' c ^ for c < 1. These states 
contain genuine three-qubit entanglement if c > (vÏ3 — 
l)/3 ps 0.869. 

Proof of Theorem 2. Based on our discussion on the 
two-qubit case, it is clear that there is a LHV model for 
von Neumann measurements on p( 3 > c ) if, and this is an 
important condition, this state is /7® 3 -invariant. After 
long but straightforward calculation, one obtains 
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(3,c) _ 



1 

8 1 



1 ® 1 



£ 5' 



— (tTfe <g> 1 ® O-Jfe 

16 



k—x,y,z 

o-k 



' CTfe ® 1) 



(12) 



This state is invariant U®U®U by inspection if we know 
that Cfc (8 <7fc is Í7 <2S £/ invariant. Thus correlation 
measurements on this state fit the LHV model given in 
Eqs. piTD]) for parties A/B, while we have (Mc)u = 
Tr(Mc\u>){u\) for C. 

We now show that p( 3 ' c ) is genuine three-partite entan- 
gled when c > (\/Ï3 — l)/3, i.e., it cannot be constructed 
by mixing pure states with two-qubit entanglement. In 
particular, it cannot be constructed by mixing different 
bipartite Werner states of the form p^ 2 ^ ' <g) pc, pa®P^ 2 ' c \ 
etc. In what follows, we adopt the definitions of Ref. [Ï3{ : 

Ri ■■= \{Wbc-V C a-V ab ), 



R-2 



{V AB - V C a), 



(13) 
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FIG. 2: Union of three sòlid ellipses: Projection of the points 
corresponding to pure biseparable states on the n/r2 plane. 
Together with the sòlid straight lines: Set of mixed bisepara- 
ble states. (diamond) State p' 3 ' 1 '. (circles) States obtained 
from 



„(3,D 



by permuting its qubits. (triangles) Two-qubit 
Werner state of the form p*- 2,1 - 1 ® 1/2 and states obtained 
from it by permuting the qubits. For more details see text. 



where \%i exchanges two qubits, and we use the nota- 
tion r/c := (Rk)- In order to examine the entanglement 
properties of the three-qubit states /o' 3 ' c \ we consider the 
projection of the set of biseparable states on the r\/r% 
plane p,3j, as shown in Fig. [21 The union of the three 
sòlid disks corresponds to the union of biseparable pure 
states of the three possible bipartitioning. One of these 
disks has the equation [Ï3| 



V3 



-n 



2V3 



r'<-:. 



(14) 



The other two can be obtained through ±120 degree rota- 
tions around the origin. The straight sòlid lines indicate 
the boundary of the convex hull of these sets. Any bisep- 
arable mixed state corresponds to a point within this set. 
Based on Eq. IT4|) and its rotated versions one can see 
that for biseparable states 



(Ri)< 



13+1 



0.77. 



(15) 



In Fig. [2] points fulfilling Eq. (fT5| are on the left hand 
side of the dashed vertical line. For p( 3 ' c ) we have ri = 
(2 + 3c) / 6 thus the state is genuine three-party entangled 
if c > (vÏ3 — l)/3. This finishes the proof of Theorem 2. 

For the state p( 3 < c ) the reduced two-qubit states pab 
and pac are entangled if c > 2/3. This means that for 
c > 2/3 a singlet between A and B, and also between A 
and C can be distilled. Hence from many copies of our 
state any three-qubit state can be obtained with local 
operations and classical communication. 
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IV. STATES p (n ' c) AND (7 0n -INVARIANT 
SYMMETRIC EXTENSIONS OF TWO-QUBIT 
WERNER STATES 

At first sight, there is no reason to limit our discussion 
to n = 3. Direct calculation, however, shows that p^ 4,1 ) 
is not £/® 4 invariant, so the previous local model cannot 
be applied. Nevertheless, the states of Ecj. (EP can be 
used to construct symmetric extensions [ÏJ, [Ï5[ . Recali 
that a (1, n— 1) symmetric extension of a two-party state 
p is an n-party operator H p such that p — r Ti3 > i,..., n Hp 
and H p is symmetric under the permutation of parties 
2,3, ...,n. For being an extension, we also need that H p 
is positive semidefmite [3] and for a quasi-extension that 
it is positive on product states @. 

It can be seen by inspection that if p( n > c > is positive 
semidefmite then it is an extension of p^ 2,c \ Symmetric 
extensions can be constructed for a larger range of c if 
the matrix p( n ' c ) is twirled: 

p (n,c) í dUU t® n pMjj^n 

Jueu(2) 

= í Af(dw)T w ®p® (16) 

Juj£C 2 ,\uj\ = 1 

where dll denotes the Haar measure, p u is defined in 
Eq. ©, and 

Tlü := \ [t ~ \ C ^ i^^k}- (17) 

k— x,y : z 

Direct calculation shows that for Werner states with p = 
2/3, 5/9, and 1/2, (l,m) symmetric extensions can be 
obtained from Eq. (fl~6|) for m = 2, 3, and 4, respectively. 

V. GENERALIZATION TO MORE PARTIES OR 
HIGHER DIMENSION 

Now rather than fixing the quantum state from the 
very beginning, we will look for the four-qubit quan- 
tum state for which the correlations fit the local model 
of Eqs. ([9l fTÜ|l . when party C and D also get physical 
qubits. Note however that there is no a priori reason why 
a LHV model should give correlations compatible with a 
quantum state. The desired state must be U® 4 invari- 
ant thus it must be a linear combination of the 4! = 24 
permutation operators [Ï3j |. It must fit all three-qubit 
correlations of our LHV model and must be invariant 
under the permutation of qubits B, C and D. It can be 
proved that this state must have the form p' :— p^A — 
K {% Hk=x,y,z < J k®°-k®cr k ®<7 k +J2i < k Il[(T k ®a k ®ai®(Ti]}, 
where II [A] denotes the sum of all distinct permutations 
of A and K is a constant. Here p^A contains the terms 
which do not affect four-qubit correlations of the form 



(u a <8> <7b ® er c ® °~d)- Setting K — 1/128 our state gives 
the same four-qubit correlations as the LHV model for 
Ma/b — °~x and M c / D = a y . Due to the finite number 
of frec parameters therc is only one such Hcrmitian ma- 
trix with unit trace. However, for Ma/b/c/d = o~x this 
matrix fails to reproduce correlations of the LHV model 
(i.e., -1/4). Thus, there is not a matrix corresponding to 
the n-qubit version of the LHV model given in Eqs. (JÜJTÜJ) 
for n > 4 qubits. 

Finally, one can explore whether the local model 
Eqs. PÍTOf. with (M c )u> = Tr(Aí c |w) can be associ- 
ated to a three-qudit state. Surprisingly, it turns out that 
the model obtained this way is not a vàlid LHV model 
for a quantum state when d > 2. In order to sce this, let 
us consider the case d — 3. Take Ma = defined by 

{a k } = {1,0,0} and {P k } = {|1)(1|, |2)(2|, |3)(3|}, and 
M' A = |1)(1|, which is actually equal to Ma but defined 
by K} = {1,0,0} and {P' k } = {|1)(1|, |2')(2'|, |3')(3'|}, 
where |2') = a\2) + 0\S) and |3') = (3*\2) - a*\3), with 
\a\ 2 + \f3\ 2 = 1. Moreover, on the other two qudits we 
measure Mb = Mc = |2)(2|. Using the methods of 
Ref. 0, we obtain (M A ® M B ® M c ) = 13/162, while 
(M' A ®M B ® M c ) = 15/162 for a = 13 = l/y/2. Thus 
(M A <8) M B «i Mc) ^ (M' A (g) M B O M c )- A similar lack 
of selfeonsisteney can be found for d > 3. 



VI. CONCLUSIONS 

We presented a family of three-qubit states for which 
correlations for all von Neumann measurements can be 
described by a LHV model. We proved that some of these 
states are genuine three-qubit entangled and distillable, 
so three-qubit entanglement is not sufficient for a state to 
be nonlocal. We also showed that there is not a quantum 
state corresponding to our model with more parties or 
higher dimension. For the details of our calculation, see 
the Appendix of Ref. [ï(|. In the future, it would be 
interesting to extend our model to general measurements. 
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Appendix 

In the Appendix we present some details of our com- 
putations. First we show how to compute the integral 
is Eq. |H|). It can easily be integrated numerically. For 
computing the integral analytically, we can rewrite it as 
an integration on the Bloch sphere 



P 



(n,c) _ 



neR 3 ,|fi|=i 



M{dÜ)gn ®p { a 



(n-1) 



(18) 



where 



8n ■= gC 1-08 ^^)' 
Pa ■= hl + Sla). 



(19) 



Here fl is the Bloch vector and sign(íl) = 
[sign(O x ), sign(ílj / ), sign(íí z )]. For computing the 
integral in Eq. (fT8| . we can use the following useful 
expressions 



neR 3 ,|o|=i 

f if m = 1 

= < S ^k=x.y,z <J k ®o k \ím = 2 
if m = 3, 



neK 3 ,|fi|=i 



M{dQ){Qa)® m (gisign(íí) 







if m = 



5 T, k =x, y ,z °k ®a k if m = 1 
if rn = 2 



4 ^-^k—x.y,z k 

J = ^ IÍ[(Jk <E) a k <E> a t (g> a{\. 



(20) 



Kk 



Here II[j4] denotes the sum of all different operators ob- 
tained from A after permuting the qubits. Based on 
these, the state corresponding to Eq. ([TT]) for n — 3 is 
Eq. (|12|) . For n — 4 the density matrix obtained from 
Eq. (HU) is 

(4,c) _ 1 (g 1 g) 1 

16 



32 



cr fe ® n[cr fc ® i ® i] 



A;— x,y,z 

— 1 ® n[cr fe ® (Tfc <g> 1] 
48 
C 

64 ( 

c 



- Y^5Z n [ ÍJfe ® cr ' £<8,cri<8,cri ]· ( 21 ) 



LHV model Eqs. (|91 1 0[) when for parties C/Z? we 
have {M C / D )w = Tr(M c/D \uj)(u;\). For p^ 4 ^ wc have 



cr-E 65 a x (8 Cu <8 fr u 



-1/8 and (aa; ® ^ ® (Tu 



— 1/4. After twirling, p( 4,c ) becomes 



p^' c) = —1(8)1(8)1(8)1 
16 



32 



cr fe (8 n[cr fe (8 1 (8 1] 



k—x,y,z 



— 1 (8 n[cr fe ® crfe (8 1] 
48 

3c 



160 

c 

Ï60 



a k <8 o-fc (8) o-fc (8 er fe 



^n[cr fe ® (Tfe (8cr; (22) 



Z<fc 

In the second part of this Appendix we show how the 
correlations for the three-qutrit case were computed. Let 
us use the notation \u>) = J2 k \fü>k exp(i0fc)|fc'), where 
u k are non-negative, 4> k are real and |1') = |1). Hence 
u k = {uj\P' k \uj). Based on Ref. 0] we can write 



(M' A (8 Mb <8 M c ) 



M(dcü)(Lü\P 2 \Lü) 2 , (23) 



where S denotes the subset of {u> € C 2 , |w| = 1} for which 
ui < ui and u\ < U3. Now, Pi in the {|2'), |3')} basis is 



Pi = 



\a\ 2 a*j3* 
a(3 \(i\ 2 



(24) 



Substituting Eq. (J24]) into Eq. (J23J) we obtain 



(M' A <8 M B <8 Ma) = / M{dw) 



\aru2 + 



+2\a(3\ y / u 2 U3 cos((/)i - 4> 3 + 4> a + <f>p) 



«3 
2 



(25) 



where cj> a and (j>p are the phase of a and /3, respectively. 
Now we transform J s M{duS) into an integral over u k and 
(f>k and do the integration [2|. Thus we obtain 



(Mi®M fl ®M c ) = ^(|a| 4 + |/3| 4 ) + g|a/3| 2 , 



whilc 
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(M A ®M B ®M C ) = —■ 
162 



(26) 



(27) 



Z<fc 



Clearly (M^ ® M B ® M c ) ^ ( M A <8 Mb <8 M c ) if < 
lal < 1. 



This matrix is not U (8 Í7 (8 U (8 Í7 invariant. How- 
ever, if for party j4 the operators Ma = & x /y/z are 
measured then the many-body correlations still fit the 



